First Excited Band of a Spinor Bose-Einstein Condensate 
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The analytical expression of the fractional parentage coefficients for the total spin-states of a 
spinor N-boson system has been derived. Thereby an S-conserved theory for the spinor Bose- 
Einstein condensation has been proposed. A set of equations has been established to describe the 
first excited band of the condensates. Numerical solution for 23 Na has been given as an example. 

PACS numbers: 03.75. Fi, 03.65. Fd 



In recent years, accompanying the experimental 
realization of the spinor Bosc-Einstein condensation 
(BEC)0,HI1I1I3, corresponding theor ies based on the 
mean field theory have been developed |rj PI IsL 1^. liH ITU 
IT2I H^ l . The total spin S together with its Z-component 
are actually conserved in the spinor condensates. Fur- 
thermore, the interactions are in general spin-dependent. 
Therefore, the strength of the mean field g that acts on 
the macroscopically occupied quantum states might de- 
pend on S. In previous theories, the S-dependence of the 
strength g has not yet been perfectly clarified. A pri- 
mary attempt along this line was proposed in Ref. [l4[ 
where the total spin-states have been introduced and the 
S-dependence of g has been derived analytically. 

The present paper is a continuation of the Ref. [l4| . 
It is noted that the thermodynamical property of con- 
densates at low temperature depends on their low-lying 
spectra, this is an interesting topic scarcely touched. The 
first aim of the paper is to generalize the study of Ref. 
[l4T | from the ground band (GB) to excited bands. Rig- 
orous total spin-states with good quantum numbers, S 
and Sz, and with a specific permutation symmetry (see 
below) will be introduced, a set of equations is derived 
to describe, as a first step, the first excited band (FEB) 
of the spinor BEC. 

When the total spin states are adopted, how to achieve 
the matrix elements among them is crucial. The second 
aim of this paper is to develop a general tool, namely the 
fractional parentage coefficients, to facilitate the theo- 
retical treatment of spinor N-boson systems. The case of 
23 Na as an example will be studied numerically. Before 
going ahead, we need some knowledge from few-boson 
systems. 

Let the interaction between a pair of spin-one bosons 
be Uij = 6{r, r ,)(),, , O tj = g m P m + g q P q = c + 
€2^1 ■ Fj . Where P m and P q are projectors to the 
two-body spin-states with spin Sij = and 2 (spin 1 is 
irrelevant since in this case both the spin wavefunction 
and the spatial wavefunction must be antisymmetric that 
prevents a point interaction), respectively, g m and g q 
are the associated strengths, Fj is the spin operator of 
particle i, c Q = (g m + 2g q )/Z , and c 2 = {g q -g m )/i. Let 
the N bosons (with mass m) be confined by a harmonic 

trap with frequency ui. When tko and y^j are used as 
units of energy and length, respectively, the Hamiltonian 



reads 

tf^EC-Vf + rfH^t/y (1) 

i i<j 

Using the method as given in |15| . the Hamiltonians 
for N = 3 and 4 two-dimensional systems have been 
diagonalized. Both the orbital angular momentum L and 
S are good quantum numbers, therefore the eigenenergies 
and eigenstates can be denoted as Els and ^ls (there 
is actually a series of states having the same L and S, 
but only the lowest of them is concerned) . It was found 
that, the levels are grouped as shown in Fig.l. In the 
lowest (first) group all bosons remain mostly in the lowest 
harmonic oscillator (h.o.) shell, in the second group one 
of the bosons is mainly in the first excited h.o. shell, etc. 
. All the levels in the first group have L — , they would 
be degenerate if the interaction were spin-independent 
(i.e., g q — g m or c<z = ), otherwise they would split into 
a band, the ground band (GB). All levels in the second 
group have L = 1 , they would be degenerate into two 
levels if g q = g m , and each would split into a band if 
9q 7^ 9m ■ Thus, the two lowest group contain totally 
three bands. In particular, we found that the third band 
is exactly a shift as a whole of the ground band by the 
unit fku , therefore it is simply the cm. excited states 
of the ground band. We shall concentrate on the second 
band, namely the FEB. 

Since the total wave function must be totally symmet- 
ric with respect to particle permutation, it can be in gen- 
eral written as 

k,Xj 

where $sf. is the total spin-state with the good quantum 
numbers S , Sz , and it is also the j — th basis function 
of the A representation of the A^-body permutation 
group. The index k refers to the multiplicity, i.e., more 
than one set of total spin states might have the same S 
and A . |Tfjj fLSk,\j is a spatial function belonging to 
the same representation A . The summation over j is 
necessary to assure *S?ls to be totally symmetric. The 
summation over A and k is necessary because that they 
are not conserved under the interaction. 
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For a given state, from (J2J we can define the weight of 
a representation as 



W LS „ 



<f. 



l/x 



LSk,\j\J LSk,\j 



> 



(3) 



When g g = g m the permutation symmetry is pure and 
therefore Wls,\ = 1 or . When g q ^ g m , the mixing 
of A appears, but remains very weak. It was found that, 
for both N = 3 and 4, the GB is dominated by A = {N} 
, while the FEB by A = {N - 1, 1}. Even if the spin- 
dependent interaction is very strong (g q /g m = 0.6 and 
9m + 9q = 3) , the GB still has W s,{n} > 0.996 and the 
FEB has Wis,{n-i,i} > 0.998 , thus the mixing is still 
weak. Since the mixing is so slight, A is nearly a good 
quantum number. Therefore, the labels (L, S, A) can be 
used to label the states as marked in Fig.l. 
There are four points noticeable. 

(i) The zero-range interaction between two bosons 
would be zero if their spatial wave function is antisym- 
metric (because in this case the two bosons do not over- 
lap). Since the interaction is repulsive, the representa- 
tion A with more antisymmetric pairs is more advanta- 
geous to binding because the repulsion can be thereby 
reduced. 

(ii) When one and only one boson has been excited to 
the second h.o. shell, the spatial permutation symmetry 
is restricted to A = {N} and {N — 1, 1}, and this explains 
why the second group contains two bands as shown in 
Fig.l. Where the band with {N — 1,1} is lower due to 
the reason stated in point (i). Thus the FEB would in 
general have A = {N —1,1} . 

(iii) The energy would also depend on S due to the 
spin-dependent interaction. Thus, the levels with the 
same A would further split into a band. 

(iv) When g q > g m , the pair coupled to Sy = 2 is 
more repulsive than the pair coupled to . Evidently, 
the many-body state with a larger S would have more 
pairs coupled to 2. Therefore, in a band, the state with 
a larger S is higher. On the contrary, when g q < g m , 
the state with a larger S is lower as shown in Fig.l. 

To generalize the above study to the case with a 
very large N, we need the knowledge of the total spin- 
states. For the most important case of A = {N}, there 

is an unique if N — S is even, or there is none 

otherwise [lj> UJ$- Therefore both the labels j and k are 
not necessary, and a simpler notation = can 
be used. 

Let us start from a 2-boson system with (12) = 
(x(l)x(2))s , where x{i) is the single boson spin-state, 
S = or 2. For N = 3 , the requirement that N — S 
is even leads to S = 3 and 1. They were given in 
as 



^ 3} = {(X(l)x(2))ax(3)} 3 = {Cx(3)}3 (4) 



Starting from J2J and JSJ, let us assume that all 
^ of the (N — l)-boson system have been known. 
Then, taking the rule of angular momentum coupling and 
the rule of outer product into account, the $g N ^ of the 
./V-boson system can be uniquely expanded as 







{N} 



{N}rAN-i} 



S+l 



X(N)} S + b { s N} {4-i 1} x(N)} 



(6) 

where and are the fractional parentage co- 

efficients (FPC) to be determined, however they are zero 
if N - S is odd. 

From 10, a similar formula can be written for $\ N ^ 



's±i 



Therefore, $g N ^ can be expanded once more as 

+ bi N} ai N - 1} {l4 N ~ 2} x(N-l)]s-MN)}s 

+ b^^H^s-fhiN-vh-ixms 



(7) 



On the other hand, since the 4^ ^ s totally symmetric 
as required, the right hand side of J7J must be invariant 
against the interchange of the N-th and (N-l)-th bosons. 
This leads to a set of three homogeneous equations. It 
turns out that they are linearly dependent, one of them 
reads 

{N},{N-1} 



[W s (S + l,S + l)-l]4 iV '& 



J S+1 



W s (S-l,S+l)b { s N} a { s N - 1} =0 



(8) 



where 



W s (T, T') = V(2T + 1)(2T' + 1)W{ISS1; TT') (9) 

and W(1SS1; TT') is the well known Racah coefficients. 
From JHl we obtain 



JN} 

a s 



7s 



dm 



{N} 

Ts 



SV2S + 1 



V2S + 3 



-/b {N - 1} 
l°s+i 



(10) 



(5+l)\/25+T 



/4-r 1} (id 



if N — S is even, where the constant 



Is 



{N} 



is introduced to 



assure (a^^) 2 + (b s N ^) 2 — 1 , so that 4^ is normal- 
ized. Eas. (jlOf) and (flip together with the condition of 
normalization lead to a unique solution, thus we obtain 
the analytical form of FPC as 

a { s N} = [(l + (-l) Ar - 5 )(7V-,5')( > 5'+l)/(2^(25+l))] 1/2 

(12) 



^ 3} = ^ 2} x(3)}r + f{4 2} x(Z)h (5) 



b { s N} = [(1 + {-l) N - s ) S {N + S+1)/{2N{2S+1))} 1/2 

(13) 
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and 



7s 



{N} 



S(S+1)[ 



(N-S)(N + S+1) 
N(N- 1)(2S + 3)(2S- 1) J 



11/2 



(14) 



The FPC are very useful because related matrix ele- 
ments among the total spin-states can be thereby calcu- 
lated. 

Now let us study the set of total spin-states with the 
{N — 1, 1} symmetry. It has been proved that the multi- 
plicity of this set is one ^(| , therefore they can be labelled 



in a simpler way as Q l s = *d 



{N-l,l}i 
S.k 



Since the dimension 



of {N — 1, 1} representation is N- 1 , i = 1, 2, • • • , N- 1 
. Making use of the FPC, they can be expanded as 

e s = ^{*SJr 1} (<)x(0}a - 4 N} {C 1} (i)x«} s 

(15) 

if N — S is even, and 



ip s can be deduced as 



,b 12 



[ho + 9s 



6 I "Sl 



where 



4nr 2 



(N — 2)g 



a I "Si 



Airr 2 



e„u 



a"s 



[h 1 + (N~l)g b s L^] u b s = e b u b s 



h - -[-— + — + r 2 ] 
2 dr 2 r 2 



(19) 



(20) 



(21) 



(22) 



(16) 



if N — S is odd, where the notation i represents all 
the bosons except the i — th . From the rule of outer 
product, l|15|) and (|16fl will contain only the symmetries 
{N} and {AT— 1,1} . However, when N—S is even, 
(|15fl is orthogonal to the unique totally symmetric spin- 
state ■ Hence, it has the pure {N—l, 1} symmetry. 

When N — S is odd, it has been stated above that t?g^ = 
, therefore (|16(l has also the pure {N —1,1} symmetry. 
It is obvious from l|12|l and (|13|l that Q l s does not exist if 
S = or N . Thus the S of @ l s is ranged from 1 to N- 1 
, while the index i of Q l s is ranged from 1 to N . This 
set of N states are not mutually orthogonal but linearly 
dependent. They satisfy ^ 85 = 0. This arises because 
the summation over i leads to a symmetrization, a state 
with A ^ {N} will not survive after the symmetrization. 

Based on the knowledge of 8^ , we can begin to study 
the FEB (three-dimensional) of realistic cases. Let the 
normalized single-particle state of condensation be de- 
noted as (fig = 

denoted as <p b s 



/ 47r r 



*4(r) 

s(r)Y h 



while the excited state be 
(0,6) with 1 = 1. As a 



basic assumption we define the spatial wave function as 



F S .i = <fs( i ) Il 3(^) l Psti) 



Then, the total wave function is written as 



*5 = E F ^ ©!< 



(17) 



(18) 



Where the summation over i leads to a symmetrization 
as required, and the {A^ — 1, l}-representation has been 
introduced. The expansion in (18) is different from (2), 
but both are eigenstates of the total spin and symmet- 
rical with respect to particle permutation, thus both are 
applicable. 

Inserting (|18|l into the Schrodinger equation, using the 
standard variational procedure, the equations for Lp a s and 



g b s = (& s \O l3 |et 



(23) 



and 



g% = (Q%\O jk |8 S ) (24) 
The energies of the states in the FEB are given by 



4 £) = e b + (N-l)e a 



(N - % 



4tt 



dr , 



(N-l)(N-2)g% 
8tt 



,a \2\„,b |2 



00 dr 

Ski 4 



(25) 



where S can be even or odd, and is ranged from 1 to 
N — 1. 

Using the expansions (|15|l and (|16fl twice (one is for the 
N — l body system), with the FPC given by l(T^|l and [tT3"jl . 

the matrix elements among 8^ can be obtained as 



9s 



N - 2 



N 



^ (0)+ c 2 )-(i + (-ir-)|e_ta_ 



c 2 (26) 



1 



9s 



Co + c 2 - 



■(N- l)(N-2) 
,N + l + (-l) N - s 



N 



S(S + 1) - 2(AT - 1)] (27) 



In the FEB, all the members with odd N — S do 
not have their partners in the ground band, they are 
the lowest states with the given S. If the system 
is cooled from higher to zero temperature, and if no 
spin-relaxation mechanism is available, the probability 
of staying in the FEB is nearly half. Incidentally, 
when N is large, since g^ is very close to the gs — 
Co + (S(S+l)-2N)/(N(N-l)) c 2 of the ground band 14 
, thus (fig is very close to that of the ground band ii N—S 
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As an application we calculated the FEB of 23 Na 
atoms trapped in a harmonic well with fku = 100 Hz 
and with g m — 6.351 x 10 _4 \/w , g q /g m = L10 . The 
normalized radial wave functions R a (r) = u a s (r)/r and 
R b {r) = u b s (r)/r are given in Fig. 2 (a) and (b) respec- 
tively. 

Let -Eg 9 ' denote the energies of the GB (E^ is the 
ground state energy ). The upper and lower bounds of 

the FEB, i.e., and E i ] relative to E ( Q 9) are plotted 

in Fig. 3. This figure demonstrates how the FEB and GB 
split. In particular when N gets larger, the lower bound 
of FEB would cross the upper of the GB, it implies an 



overlap of these two bands. 

In summary, we have derived a set of equations to de- 
scribe the first excited band, this might lead to a deeper 
understanding of the condensates at very low temper- 
ature. Furthermore, we have achieved the analytical 
form of the fractional parentage coefficients. This will 
greatly facilitate the calculations dealing with the total 
spin-states, and will thereby promote the applicability 
of the S-conserved theory for various spinor many-body 
systems. 
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FIG. 1: The evolution of Els (in hu> ) with g q /g m in 3 
-boson (a) and 4-boson (b) systems. (L, S, A) are marked by 
the levels. 
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FIG. 2: The radial wave functions, R a (in (a) ) and R b ( in 
(b) ), as functions of r. The dashed curves are for N = 100, 
while the solid ones are for N = 10 5 . In each case the curves 
for S = N — 1, N — 2, 2, and 1 are plotted. However, they 
can not be distinguished except the case of R b with a large 
N, where S are marked by the curves. In (b), when S lies in 
between and N — S is even (odd), the associated curves are 
distributed between the curves with S = 2 and 99998 (1 and 
99999 ). Obviously, when N — S is odd, the S-dependence is 
less explicit, and the R b with S = N — 2 has a smallest size. 




FIG. 3: E { s e > - E { 9 ' against N. The upper (lower) solid line 
has S = N - 1 (1). The dashed line is for - E^ ] , namely 
the split of the GB. 



